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Abstract 

First correction to the high-energy asymptotics of the total e + e~ photoproduction cross section 
in the electric field of a heavy atom is derived with the exact account of this field. The consideration 
is based on the use of the quasiclassical electron Green function in an external electric field. The 
next-to-leading correction to the cross section is discussed. The influence of screening on the 
Coulomb corrections is examined in the leading approximation. It turns out that the high-energy 
asymptotics of the corresponding correction is independent of the photon energy. In the region 
where both produced particles are relativistic, the corrections to the high-energy asymptotics of 
the electron (positron) spectrum are derived. Our results for the total cross section are in good 
agreement with experimental data for photon energies down to a few MeV . In addition, the 
corrections to the bremsstrahlung spectrum are obtained from the corresponding results for pair 
production. 
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I. INTRODUCTION 



Knowledge of the photoabsorption cross sections is very important in various applications, 
see, e.g., jl|. The relevant processes are the atomic photoeffect, nuclear photoabsorption, 
incoherent and coherent photon scattering and e + e~ pair production. In the coherent pro- 
cesses, by definition, there is no excitation or ionization of an atom. The high-accuracy 
estimation of the corresponding cross sections is required. They have different dependence 
on the photon energy uj. At u > lOMeV, the cross section of e + e~ pair production becomes 
dominant Q] • The coherent contribution a co h to the pair production cross section is roughly 
Z times larger than the incoherent one [Z is the atomic number), thereby being the most 
important for heavy atoms. Just the coherent pair production is considered below. 

The theoretical and experimental investigation of the coherent pair production has a 



long history, see 



2 



_ . In the Born approximation, the cross section as is known for arbitrary 
photon energy flfl. The acconnt of the effect of screening is straightforward in this ap- 
proximation and can be easily performed if the atomic form factor is known jf|. For heavy 
atoms it is necessary to take into account the Coulomb corrections ac, 

O'coh = o B + o c ■ (1) 

These corrections are higher order terms of the perturbation theory with respect to the 
atomic field. The magnitude of ac depends on uj and the parameter Za (a = 1/137 is the 
fine-structure constant). The formal expression for ac, exact in Za and uj, was derived by 
0verb0 et al. jf]]. This expression has a very complicated form causing severe difficulties 
in computations. The difficulties grow as uj increases, so that numerical results in 6] were 
obtained only for uj < 5MeV. 

In the high-energy region to ^> m (m is the electron mass, h = c = 1), the consideration 
is greatly simplified. As a result, a rather simple form was obtained in [3,0] for the Coulomb 
corrections in the leading approximation with respect to m/u. However, the theoretical de- 
scription of the Coulomb corrections at intermediate photon energies (5 4- lOOMeV) has not 
been completed. At present, all estimates of ac in this region are based on the "bridging" 
expression derived by 0verb0 [jj. This expression is actually an extrapolation of the results 
obtained for uj < 5MeV. It is based on some assumptions on the form of the asymptotic 
expansion of ac at high photon energy. It is commonly believed that the "bridging" ex- 
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pression has an accuracy providing the maximum error in a co h of the order of a few tens of 
percent. 

Here we develop a description of e + e~-pair production at intermediate photon energies 
by deriving the next-to-leading term of the high-energy expansion of ac- First we consider 
a pure Coulomb field and represent <Jc in the form 

a c = of + aff + 4 2) + • • • (2) 



The term Oq has the form {m/u) n S^(\auj/m), where S^ n \x) is some polynomial. The 



UJ- 



independent term affl corresponds to the result of Davies et al. |8(. In the present paper we 
r (1) 
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derive the term Oq. It turns out that is tu-independent in contrast to a second-degree 
polynomial suggested by 0verb0 j^. We present an ansatz for which provides a good 



agreement with available experimental data for uj > 5MeV. 

The high-energy expansion of the Coulomb corrections to the spectrum has the form 
similar to (J2J). In the region e± ^> m, we derive the term da^ /dx, where e_ and e + are 
the electron and positron energy, respectively, x = e_/uj. The term da^/dx may turn 
important, e.g., for description of the development of electromagnetic showers in a medium. 
The correction found is antisymmetric with respect to the permutation e + «-> e_ and does 
not contribute to the total cross section. In fact, originates from two energy regions 
e + ~ m and ~ m, where the spectrum is not known. However, our result for cr^ allows 
us to claim that the spectrum in these regions differs drastically from the result obtained 
by Davies et al. [8| for e± 3> m, if the latter is formally applied at e_ ~ m or e + ~ m. 

The effect of screening on <jq at uj ^> m is considered quantitatively. In the leading 
approximation, we find the corresponding correction (Tq CV \ which is a;- independent similar 
to a c . So, for the atomic field, a c should be added to the right-hand side of Eq. (JzJ). 
The screening correction to the spectrum is also obtained. 

In this paper we present the explicit calculations of the corrections, which have been 

n 

given without derivation in our recent work 10] and used for the detailed comparison of 
theory with experimental data. 
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II. GENERAL DISCUSSION 



The cross section of e + e pair production by a photon in an external field reads 

d(T coh = dpdq5(u - e + - e_)\M\ 2 , (3) 

where e_ = e p = \/p 2 + m 2 , e + = e q , and p, q are the electron and positron momenta, 
respectively. The matrix element M has the form 

M = J dr^ p +) (r) e tp^ (r) exp (ikr) . (4) 

Here i/j^ and ip q ^ are positive-energy and negative-energy solutions of the Dirac equation 
in the external field, e M is the photon polarization 4-vector, k is the photon momentum, 
e = e^j^, 7 M are the Dirac matrices. It is convenient to study various processes in external 
fields using the Green function G(r 2 ,ri\e) of the Dirac equation in this field. This Green 
function can be represented in the form 

G(r 2 , ri|e) = 



4 + y 2 )4 +) (n; 

e — e n + iO 



+ 



dp 



e — e p + iO e + Bp — iO 



(5) 



where ipn is the discrete-spectrum wave function, e n is the corresponding binding energy. 
The regularization of denominators in (jSJ) corresponds to the Feynman rule. From (jHJ), 

.(2n) 2 



\2 



-5G (r 2 , ri\ - e 



I dO p ^ +) (n)^ +) (r 2 ) = i^SG( ri ,r 2 \e p ), (6) 

where Q p is the solid angle of p, and 5G = G — G. The function G is obtained from © by 
the replacement iO <-> — iO. 

Taking the integrals over Q p and Q q in (JHJ), we obtain the electron spectrum, which is 
the cross section differential with respect to the electron energy e_ . Using relations (jSJ, we 
express this spectrum via the Green functions: 

(]<J "' 1 ' -!! d ri dr 2 e- lkr Sp {SGfanleJeSGin,^ ~e + )e} , (7) 



de- u> _ 

where r = r 2 — V\ and e + = uo — is the positron energy. Since the spectrum is independent 
of the photon polarization, here and below we assume e* = e (linear polarization). 



Due to the optical theorem, the process of pair production is related to the process of 
Delbruck scattering (coherent scattering of a photon in the electric field of an atom via 
virtual electron-positron pairs). At zero scattering angle, the amplitude Mp of Delbruck 
scattering reads 



M, 



D 



2ia Jde J J dr 1 dr 2 e ikv Sp {G(r 2 , r^e) eG(r x , r 2 \e - u) e] . (8) 



It is necessary to subtract, from the integrand in (JBJ), the value of this integrand at zero 
external field (Za = 0). Below, such a subtraction is assumed to be made. 

It follows from Eqs. ©,© and the analytical properties of the Green function that 

— Im M D = a coh + a bf . (9) 

Here 

o- 6/ = — % 1^L if dr 1 dr 2 e~ tkr ^2 Sp{p n (r 2 ,r 1 )e5G(r 1 ,r 2 \e n - u)e} , 

Pn(r 2 , ri) = lim (e - e n )G(r 2 , n\e). (10) 

The quantity a b f coincides with the total cross section of the so-called bound-free pair 
production when an electron is produced in a bound state. In fact, due to the Pauli principle, 
there is no bound- free pair production on neutral atoms. Nevertheless, the term a b f should 
be kept in the r.h.s. of Q. In a Coulomb field, the total cross section a b f was obtained in 



111 ] for u ^> m. In this limit, a b f oc 1/mu and should be taken into account when using the 
relation @ for the calculation of the corrections to a co h from. The main contribution to <j b f 



comes from the low-lying bound states [ll| when screening can be neglected. So, in we 
can use a b f obtained in 

It is convenient to represent do cobi jde^ and Md in another form using the Green function 
D(r 2 ,ri\e) of the squared Dirac equation, 

G(r 2 ,n|e)= [ 7 (^-V r (r 2 ))-7p 2 + m] J D(r 2 ,r 1 |£), p 2 = -iV 2 (11) 

n 

According to |l2j, we can rewrite Eq. (J7J) in the form 
_ « ff dri dr 2 e-^ 



ds- 2uj 

xSp{[(2ep 2 - ek)6D{r 2 , n|e_)][(2epi + ek)5D(r u r 2 \ - e + )]} , (12) 
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and Eq. (JHJ) as 

M D = ia de dr 1 dr 2 e~ ikr 



xSp{[(2ep 2 - ek)D(r 2 , n\u - e)][(2e Pl + ek)D(r u r 2 \ - e)}} 

+2ia de dr SpD(r,r\e) . (13) 



The last term in (J13)) is ^-independent, and has no imaginary part. Therefore, it does 
not contribute to the relation (jJJJ). 

III. GREEN FUNCTION 

To obtain the spectrum ((Zj), (|T2^1 and the Delbriick scattering amplitude (JBJ), (JEJ) it is 

necessary to know the explicit form of the Green function of the Dirac equation in the 
Coulomb potential V(r) = —Za/r. An integral representation for G(t2,ti\e) has been 
obtained in For |e| > m it has the form 

i f°° 

G(r 2 ,ri\e) = / dsexp[2iZas A + w(r 2 + ri) coth s] T 

\~r>f\K Jo 

T = [1 - (7 • n 2 )( 7 ■ n 1 )][(7°e + m)| d y S B - iZa^KCoths S B ] 
+ [1 + (7 ■ n 2 )( 7 • n x )](7 £ + m)S A + imZa^-f ■ (n 2 + uJSb 

-7 • (n 2 +7ix)S B — /«coths 7 • (n 2 - ni)^. (14) 



2 sinh 2 s 



In this formula 

00 00 
5a = ^e-^J 2 ,(y) Z[P/(x) + P/^rr)] , 5 fl = ^e^ J 2 ,(y)[P/(x) - PU(x)} 



1=1 

v = \Jl 2 — (Za) 2 , k = \Je 2 — m 2 , X — e/k 



y = 2K v /r^r 2 "/sinhs, x = n x • n 2 , m, 2 = r lj2 /r x>2 , (15) 

J2u(y) are Bessel functions and P(x) are Legendre polynomials, P((x) = d x Pi(x). The Green 
function D(r 2 , rx\e) can be obtained from ((TH) by keeping in T the terms oc m: 

i f°° 

D(r 2 ,ri\e) = / dsexp[2iZas A + iK{r 2 + ri) coths] 

Anr 2 r 1 K J 

x { [1 - (7 ■ n 2 )( 7 • nx)]| d y S B + [1 + (7 ■ n 2 )( 7 ■ nOJSU 
+zZa7 7-(n 2 + n 1 )5 , B }. (16) 



We are going to derive the high-energy asymptotic expansion of the spectrum in the 
region e± ^> m. For the first two terms of such an expansion the main contribution to 
the integral in (J7J), (fT2|) is given by the region r = \r 2 — fi| ~ u/m 2 , see 3]- Let us 
introduce the variable p as the component of ri (or r 2 ) perpendicular to r: 

p= rxj£xr>] (i7) 

As shown in Q|, the main contribution to the Coulomb corrections to the spectrum 
originates from the region p ~ 1/m and 9, ip ~ m/cj <C 1, where 9 is the angle between the 
vectors r 2 and — r\, and ?/> is the angle between the vectors r and k. In this region we have 
9 m rpjr\r 2 . The argument of the Legendre polynomials in ()15|) is x = rii ■ n 2 ~ — 1 + 9 2 /2. 
Besides, the term k(t 2 + ri) ~ uj 2 /m 2 ^> 1 in the exponents in (|14|). (|T^|) is large, and 
the integral is determined by large s. Then coths ~ 1 + 2exp(— 2s), and from (|T6|) we 
have exp(— 2s) ~ 1/rer ~ m 2 /uj 2 . The argument, y, of the Bessel functions in Sa,b can be 
estimated as y ~ kt/ sinhs ~ tu/m ^> 1. 

A simple method of the calculation of Sa,b at y 3> 1, 1 + a? ~ # 2 /2 <C 1, and y9 ~ \ has 
been formulated in the Appendix of ^|. It turns out that the leading term and the first 
correction are determined by values of / ~ y ~ uo/m in sums Sa,b ■ This fact is in agreement 
with the evident estimate I ~ ep ~ uo/m ^> 1. In the same way as in we obtain for Sa,b 
with the first correction taken into account 



MyO/2) 



l + i 



S 



B 



y_ 

29 



MvB/2) 



1 + 



7T 



(Za) 



;is) 



Let us pass in (fTH|) from the integration over s to the integration over y, see (JT5J). We have 



[2iZaXs] n ( W^ V^ coths 



cxp 



\ y J 



i + 



\k 2 T\T 2 



Then we obtain 

# fa, file) 



i/t(ri+r 2 ) 

\^-nKr 2 r\ j 



ydy 



4/t A /r 1 r 2 



2iZaA 



cxp 



x|[l+m(Za) 2 /y] 
+7r(Za) 



8/trir2 

J (j/g/2)+2iZ Q! ' 7l( y ) a-(n 2 



yO 



2 



n 2 x ni] • £ j , 



(19) 



where en = 7°7, £ = (z/2)[7 x 7]. This expression is the quasiclassical Green function of 
the squared Dirac equation with the first correction taken into account. The leading term 



in this expression, as well as the corresponding expression for G(r 2 ,ri\e) has been derived 



in 



14. 



15|. It is convenient to rewrite (j!9|) in another form using the relations 



qdqJ (q9)g(q)= / ^ e g(q) 



dq 



qdqJ 1 {q6)g(q) 



Uq(q-0) ■ e 



9(q), (20) 



2tt qd 

where g(q) is an arbitrary function, q and 6 are two-dimensional vectors. In our case we 
direct the vector along p so that 6 = rpjr\r 2 . Using (J2(J|) we have 



D(r 2 ,rx\e) 



^gi/t(ri+r- 2 ) 



57TKr 2 ri 



dq 



q 



2iZaX 



exp 



1 + 



2g 



Zoi 



ct ■ q 

q 2 



i[r 2 + ri)<r 
2nr 1 r 2 
in(Za) 2 



+ iqd 



r/r)-[qxE]}. (21) 



4g 3 

The leading term of this formula has been obtained in ^|. Let us integrate by parts 
the term containing o;-matrix and make the change of variable q — > n(q — p). We obtain 



■D(r 2 ,ri|e 
1 + 



iKe 



i>7r*rir 2 



x 



2r\r 2 X 



(X ■ q 



dqexp 
l + i 



nrq 
2r x r 2 
ir(Za) 2 



2n\q 



2jr^F 2 

\Q ~ P\ 

7r(Za 



2iZa\ 



(22) 



An 2 



.( 7 o/ A _ 7 . r/r ) 



7 • (q - p) 



Note that in this formula and below we can set A = signs. 

It is easy to check, that within our accuracy the contribution of the last term in braces 
vanishes after taking the trace in (|12p. Therefore, this term can be omitted in the problem 
under consideration. The remaining terms in (J22j) can be represented in the form 

at - (pi + p 2 ) 



D{r 2 ,r x \e) = 
D^{r 2 ,r x \e) 



1 + 



2s 



D(°\r 2 , ri \e) 



me 



87r 2 rir 2 



dqexp 



nrq 
2r x r 2 



2^r~ 2 



2iZaX 



l + i 



7r(Za) 
2n\q — 



(23) 
(24) 



The function D^(r 2 ,r 1 \e) is nothing but the quasiclassical Green function of the Klein- 
Gordon equation in the Coulomb field. The function 5D in (|T2*|) is defined as SD = D — D, 
where D is obtained from ()23|) by the replacement — > D^*. 



IV. COULOMB CORRECTIONS TO THE SPECTRUM 

In this section we consider the Coulomb corrections to the spectrum, dac/dx, for e± ^> m 
taking into account terms of the order m/e±. According to |8[, the higher order terms of 
the perturbation theory with respect to the external field (Coulomb corrections) are not 
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seriously modified by screening. However, this question has not been studied quantitatively 
so far. The influence of screening on Coulomb corrections is investigated in detail in Section 
IV1I In the present Section we calculate doc / de_ in a pure Coulomb field. 
Substituting (J23j) in (J12j) and taking the trace, we obtain 
dac 4a , 



-Re J J dr 1 dr 2 e- lkr ^A[e-p 2 D W ][e-p 1 D 



f] 



D 



(°) ^ D (o) 



[e-(p 1 +p 2 )D W ][e.( Pl +p 2 )D^]} 



(0) 



(r 2 ,ri|e_j, u+ = n^{r u r 2 \ - e+) . (25) 

The terms oc D^D^* are omitted in this formula since they do not contribute to the leading 
term and the correction we are interested in. Besides, we have integrated by parts the terms 
containing second derivatives of D^°\ In this formula and below we assume the subtraction 
from the integrand the terms of the order {Za)° and (Za) 2 . Then we use the relation 



in 2 e tKr 



(e-p 1>2 )D^(r 2 , ri \s) = ^~^rj dqexp 

x 1 + z — } '— 

\ 2n\q — p\ 

and pass from the variables r\, 2 to the variables 

r x [t*i x r 2 ] 
r = r 2 -r h p = - , 



nrq 2 1 
2rir 2 
e ■ r 



T- 



r 



\Q ~ P\ 

r l,2 



2iZa\ 



(26) 



K r ■ r l, 



(27) 



In terms of these variables dr\dr 2 = r dr dp dz, and within our accuracy r\ = rz, r 2 
r(l — z). We obtain from (J25|) 



da, 



c 



cfe_ 
x exp 

x <^ Ass 



ae_e + 



Re 



dr 



dz 



luir 



m 



e ■ r 



E-E + 

e • q 



1 



z 



z 2 (l-z) 
e-ql+e + q 2 + 
% 2rz(l-z) 
e ■ q 



dq^dq + dp 

in(Za) 2 



0, 

Q- 



i + 



2iZa 



1 1 

+ 



-e ■ r + 



£-Q- e+Q 

: , {] _ z]2 ,e-q-)(e.q + ) 



uj- 



(28) 



where Q± = \q± — p\ and ip is the angle between vectors r and k. Since dac/de- is 
independent of the photon polarization, we can replace in (f2~H|) e % e? by \5 l [ = |(<5* J — k 1 ^ /uj 2 ). 
The integral over p can be taken with the help of the relations (see Appendix) 

2iZa 



f(Za) 



g(Za) 



2ir(Za) 2 q 2 



4nq J Q 



dp 



dp 

Q- 



Q± 

Q- 

2iZa 



-l + 2(Za) 2 ln 2 



0, 

Q- 



Za 



= Re[pl + iZ a) + C] 
r(l - iZa)T(l/2 + iZa) 



r(l +iZa)T(l/2-iZa) 



(29) 
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where ip(t) = d In T (t) /dt, C = 0.577... is the Euler constant, q = |q_ — q + \. We have 
da c a(Za) 2 e-e + f°°dr f°° , f 1 dz 

WJJ aq ~ 

g(Za) g*(Za) 



de- 8ti 2 uj J r 3 J J z 2 {l — z 



x exp 



1 , Re / ^3 / ^ / y 2M _ ^2 // d 1- d( l 
i?- — J I 



iur ( 2 m 2 \ .£_g 2 + £+g ,2 _ 



2 V £-£+/ 2rz(l - z) 



q f(Zai) + 7rg 



- " lc - c+ ( ~ rV + Jr^y ) _ ^ ( r_^2 (g- - <?-!-> !- ■ (30) 



,2 



Passing to the variables q = q + q + , q = q — q + , we take all integrals in the following 
order: dip, dq, dq, dr, dz. The final result for the Coulomb corrections to the spectrum 
reads 



da$ da^ 



+ —r- = -4a 



dx dx 



l-^x(l-x))f(Zn) 



vr 3 (l - 2x)m A3 , , 

1 1 — -x[l — x) Reg (Zee) 



8x(l-x)u V 2 
x = e^/u , cr = a(Za) 2 /m 2 . (31) 

In (|31j). the term oc f(Za) corresponds to the leading approximation da^ /dx j^J, the term 
oc Reg(Za) is the first correction da^ /dx. In contrast to the leading term, this correction 
is antisymmetric with respect to the permutation e + <-> (or x «-> 1 — x) and, therefore, 
does not contribute to the total cross section. Besides, the correction is an odd function 
of Za due to the charge-parity conservation and the antisymmetry mentioned above. The 
antisymmetric contribution enhances the production of electrons at x < 1/2 and suppresses 
it at x > 1/2. Evidently, the opposite situation occurs for positrons. Qualitatively, such a 
behavior of the spectrum takes place for any u being the most pronounced at low photon 
energy At intermediate photon energies, the spectrum (|3~Tj) essentially differs from that 
given by the leading approximation. We illustrate this statement in Fig.^ where a^ l dac / dx 
with correction (solid line) and without correction (dashed line) are plotted for Z = 82 and 
u = 50 MeV. 

Due to the antisymmetry of da^ /dx at e± ^> m, the term cr^ in the total cross section 
may originate only from the energy regions ~ m and e + = u — E- ~ m. The quasiclassical 
approximation can not be used directly in these regions, and another approach is needed 
to calculate the spectrum. We are going to do this elsewhere. However, for the total cross 
section, it is possible to overcome this difficulty by means of dispersion relations (see Section 
El). 
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FIG. 1: The dependence of a Q dac/dx on x, see (|31[). for Z = 82, a; = 50MeV. Dashed curve: 
leading approximation; solid curve: first correction is taken into account. 

As known [see, e.g.,[l8f, the spectrum of bremsstrahlung can be obtained from the spec- 
trum of pair production. This can be performed by means of the substitution e + — > — e 
, oj — > —oj'i and dx — > ydy, where y = o//e, u' is the energy of an emitted photon, e 
is the initial electron energy. Using (|31|). we obtain for the Coulomb corrections to the 
bremsstrahlung spectrum 



V 



dal 
dy 



-4<7n 



7T 3 (2 — y)m 



(32) 



8(1 -y)e 

This formula describes bremsstrahlung from electrons. For the spectrum of photons 



emitted by positrons, it is neces sary to change the sign of Za in (|3*2*|) . Our result (|3*2*|) 
coincides with that obtained in if the obvious mistake in the latter is corrected by 
changing 



1 1 

7^2 

The correction 



m m 
— + 

e e 



to' 



(2 - y)m 
2(1 -y)e 



in Eq.(22) of Q|. The correction (|3*2*j) is the most important at y close to unity, see Fig. 
I2J where (Tq yda^/dy with correction (solid line) and without correction (dashed line) are 
shown for Z = 82 and e = 50 MeV. 
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FIG. 2: The dependence of GQ 1 yda 1 c /dy on y, see (|3Tj) . for Z = 82, e = 50Mey. Dashed curve: 
leading approximation; solid curve: first correction is taken into account. 

V. COULOMB CORRECTIONS TO THE TOTAL CROSS SECTION 

In the leading approximation, the Coulomb corrections, to the total cross section 
of pair production for uj ^> m were obtained in j^j]. Using this result and dispersion rela- 
tions, the corresponding term, M<°> , in the forward Delbrnck scattering amplitude M D was 
obtained in [20]. These two quantities read: 

4 0) = -f ^ f{Za) , M£l = -i^uao f(Za) , (33) 

where <To and f(Za) are defined in (}3"Tj) and (|2*9*jl . respectively. 

In this Section we derive the correction a$ by means of the relation ([9*]). Starting from 
(fT3j) and performing the same calculations as in the previous Section we obtain 

M§1 + M^ c = -Aiuoa jf dx ^1 - ~z(l - x)^ /(2T«) 

_^( 1 _| I(1 _,)(^_^)) ] . 

Here the integration over x corresponds to the integration over e/uj. After the integration the 
in (|34|) coincides with that in ()33|) . The integral in is logarithmically divergent. 
Note that we have obtained the integrand in (|3*3*|) under the conditions x ^> m/o; and 
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1 — x ^> m/uj. Taking the integral from S to 1—5, where 5 > m/u, we find within 
logarithmic accuracy that Im M^ c vanishes and 

Re M%1 = <*V<*?****9{Z«) ln u 

DC m m y ' 

The quantity ImM^. does not contain ln(u;/m) and is determined by the regions of inte- 
gration over e, where e ~ m and uj — e ~ m, and, therefore, the quasiclassical approximation 
is invalid. Nevertheless, this quantity, which is related to (Jq ©, can be obtained from the 
dispersion relation for Mr, |20] 

t~» n jt / \ 2 2 „ Z" 00 ImM D (u/) duj' 

ReM D (u) = -uj 2 P DK ; . 36 

7T UJ\UJ U — UJ 1 ) 

Using this relation, it can be easily checked that the high-energy asymptotics (|35|) unam- 
biguously corresponds to the ^-independent high-energy asymptotics 

ImMS = - aN ^ mgN ' (37) 
Substituting (fSTj) into © and using 05/ from in the form 

a bf = Arra (Za) 3 h{Za) - , (38) 

UJ 



(i) 



we have for a c 



' 4 

4 1} = -*o yIm^(Za)+4vr(Z«) 3 /i(^) 



m 



(39) 



The function f\(Za) is plotted in Fig. El 

The quantity {u/m)a c </oq is shown in Fig. 0] (solid curve). It is seen that this ratio 
is numerically large for any Z. Therefore, the term <Tq gives a significant contribution to 
ac for intermediate photon energies. Dashed curve in Fig. H] gives the same ratio when cr^ 
in (J3~9"j) is omitted. It is seen that the relative contribution of the term oc f\{Za) in (J3T?|) is 
numerically small. 



VI. SCREENING CORRECTIONS 

In two previous Sections the cross section of e + e~ pair production has been considered for 
a pure Coulomb field. The difference, 5V(r), between an atomic potential and a Coulomb 
potential of a nucleus leads to the modification of this cross section known as the effect of 
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FIG. 4: The quantity (u/m)(TQ /(Tq as a function of Z(solid curve). Dashed curve corresponds 
to the same quantity without the contribution of the bound- free pair production. 



screening. In the Born approximation, this effect was studied long ago [see, e.g.,[5(. Let 
us consider now characterizing the influence of screening on the Coulomb corrections. 
Recollect that the Coulomb corrections denote the higher-order terms of the perturbation 
theory with respect to the atomic field. So far it was only known that the correction 
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is not large [8|. Here we consider this issue quantitatively. 

The quasiclassical Green function D^ ' (rz, ri\e) for an arbitrary localized potential V(r) 



has been obtained in 



the form (see also 



21] with the first correction taken into account. The leading term has 



D( \r 2j rd e) = — / dqexp i iXr / dxV (ri + xr — q 

87r 2 rir 2 J I 2r x r 2 J 

Substituting this formula into (J23J), we obtain (cf. (12*%!) ) 

p. = _^£±Re f % f jjj d M p 

cte_ 16vr 4 u; J r 5 Jo z 2 {l-z) 2 JJJ 



(40) 



i 2 (l -z)< 



(41) 



x {4^+ (e • r + ^) (-e • r + - — (e • 9 _)(e ■ ^ 

$ = r / dx[V(ri + xr — qr + ) — V(ri + xr — q_ 
Jo 

The phase $ can be represented as 

$ = 2Zaln(Q+/Q_) + $( scr ) 

= 2Za m(Q + /Q_) + r / dx[<JV(ri + xr - q+) - (JV(ri + xr - q_)\ . (42) 

Jo 

As in the case of a pure Coulomb field, the main contribution to the Coulomb corrections 
comes from the region of integration q± ~ p ~ 1/m. The main contribution to the integral 
over x in (}4*2*|) comes from the narrow region around the point Xq = —r\ ■ r /r 2 = z, 5x = 
p/r <C 1. Therefore, it is possible to perform the integration in (}4*2~j) from — oo to oo. Thus 
we can estimate <3>( scr ) as $( S£T ) ~ p5V(p) ~ Za<5V(p)/V(p) <C 1. In our calculation of 

T (sc 



(?Q Cr \ we retain the linear term of expansion in $( scr ). By definition, 



where F(Q) is the atomic electron form factor. Substituting this formula to (|42|) and taking 
the integral over x from — oo to oo, we obtain for <3>( scr ) 

$(-) = / ( e «A ± (p- g+ ) _ e iA x(p - q _)) F ( A± )^ , (44) 

where is two-dimensional vector lying in the plane perpendicular to r. Then we use the 
identity, see Eqs. (22), (23) in |l7| 

2iZa 



//I I \ AlZjOL 

dp ex P l iAl - ' (P ~ 9±)1 
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4A^ 



df U±j exp [iq ■ f ± /2] 



(45) 



where q = q — q + and f± = f ± Aj_. Expanding (JIT]) in $( scr ) and using the identity (|43j). 
we take the integrals over variables q±, r, and z and obtain 



da, 



(scr) 
C 



dx 
R(fi, a) 



8a(Za) 
3tt 

(m-i; 



1 



R(£-,a) R(£ + ,a) 



201 



4/x 2 

—18 — a(/i — 3) 
£± = 1 + 16m 2 //l , a = 6z(l - x 



2vr V/- 
[18 - 6/i + a(/i 2 + 2/i - 3)] In 



f 2 - _ ft 



.01- 1. 



(46) 



Using the trick introduced in Section IX in [17| , we rewrite this formula in another form. 
Let us multiply the integrand in (J4T)|) by 



dyS (^y 2 f A 



(/ 2 + A 2 , 



1 dy 



-l \y\ 



S((f - A ± / y y - Ai(W - i)) , 



(47) 



change the order of integration over / and y, and make the shift / — > f + A.±/y. After that 
the integration over / can be done easily Then we make the substitution y = tanhr and 




obtain 



da. 



(scr) 
C 



32 



dx 3 

r2 * dip 



dQ 

Q 



3 F(Q) 



dr 



sinh r 



sm(2Zar) 
2Z^ 



2tt 
1 + 



e~ T R{n- 



lm 2 e ±T sinh 2 r 



Q 2 (coshr + cos<^) 
Integrating over x, we have 

( SC r) 32 2 Z" 00 r 

jf " ^ [e r i?(// + , 1) - e-^(/i. 



(48) 



dr 



c 



sinh t 



sin(2Zar) 
2Zo< 



— T 



(49) 



Similar to (Tq , this correction is ^-independent. Shown in Fig. is the Z- dependence of the 



'c 

ratio a^^/a^ calculated with the use of the form factors taken from |22jJ. As seen from 
Fig. this ratio is approximately fitted by the linear function, a 



(scr) 
C 



-5.4 ■ 10~ 4 ■ ZaP. 



The corresponding correction to the bremsstrahlung spectrum is obtained from ()48j) by 
means of the same substitutions as in Section ITVl So that the quantity y~ 1 da^ sc ^ jdy is 
given by the right-hand side of (jlEJ) if we set a = 6(y — l)/y 2 - 



VII. ESTIMATION OF aff FROM EXPERIMENTAL DATA 

The most detailed and accurate experimental data have been obtained just in the region 
of intermediate photon energies. In this region, the first correction <Jq\ obtained above, 
becomes large, see Fig. H]for (uj/m)cr^ /(Tq\ and the next term Oq in the expansion (J2J) 
may be significant. Using the arguments similar to those presented by Davies et al. [8] the 
following ansatz for cr^? has been suggested in our recent paper 

aP = a [b\n(u/2m) + c] (^j* , (50) 

where b and c are some functions of Za. It was shown in that experimental data for 
<y co h are well described by the formula 

(0) , (scr) (1) (2) / K1 s 

Vcoh = + (?c + °C + a C + a C i ( 51 ) 

where Gq\ <J^ cr \ and are from Eqs. (f3~3*)). (jjHj), and (jHHJ), respectively; is given by 
flUJ) with b = -3.78(a;/m)cr ( 7Vg ) , c = 0. 
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It is interesting to compare ourpredictions for the Coulomb corrections to the total cross 



section with the results of 0verb0 



Shown in Figs. |H]and[7|is the ratio S 
r (o) 



{&coh-(7B)/cr, 



(0) 
C i 



which is the Coulomb corrections in units of Our results are represented by 





the experimental data are also shown. The results for Bi are plotted in Fig. 13 with the 



experimental data taken from 



experim ental dat a taken from 



those of 



23 L The results for Pb are plotted in Fig. [7| with the 
2J, |25( . It is seen that the difference between our results and 



)verbd is small at relatively low energies and becomes noticeable as u increases. 



According to our results, this difference tends to a constant Cq /uq^ at 



10 



oo. The 



ex perimen tal data are, on the whole, in a better agreement with our results than with those 



of 



0verbo. 



VIII. CONCLUSION 

For the e + e~ photoproduction, we have calculated the leading correction 1)31)1 to the 

electron spectrum in the region e± ^> m. This contribution noticeably modifies the spectrum 

at intermediate photon energy. It turns out that the correction is antisymmetric with respect 

to the permutation e + <-> e~ and hence does not contribute to the total cross section. The 

leading correction to the total cross section, (Tq , originates from two regions, e + ~ m and 

We have obtained q® flU using dispersion relations. In contrast to the form of 

the fit suggested by 0verb0 [9J, the quantity Oq does not contain any powers of ln(u;/m). 

We have also performed the quantitative investigation of the influence of screening on the 

Coulomb corrections ()48j) . (J49)) . It is important that does not vanish in the high- 

ly) 

energy limit. We have suggested a form for the next-to-leading correction, a c , to the total 
cross section. Altogether, the corrections found allow one to represent well the available 
experimental data. 

Starting with the results obtained for the e + e" photoproduction spectrum, we have ob- 
tained the corresponding corrections to the bremsstrahlung spectrum as well. 
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APPENDIX 



In this appendix we derive the formulas (|29|) . In the integral for f(Za) let us make the 
change of variables p — > (p + q + + q~)/2: 

2iZa 



f(Za) 



1 



87r(Za)V 



dp 



|P + g| 

\p-q\ 



l + 2(Za) 2 \n 2 lP+ql 



\P-Q\ 



(Al) 



Let us multiply the integrand in (fltjj) by 



dySly 



-1 



2pg 

p 2 + g 2 



p 2 + <z 2 ) f p-M P -q/y) 2 - q 2 (i/y 2 -i)) 
J-i \y\ 



(A2) 



change the order of integration over p and y, and make the shift p — > p + q/y. Then the 
integral over p becomes trivial and we have 



f(Za) 



1 



A{Zaf J_ x ]yj 



i-y/ 2 Vi -i/ 



2/1 + 2/ 



(A3) 



Then we make the substitution y = tanh r and obtain 

1 



f(Za) 



2{Za) 2 



o 



(ire 



sinh t 



00 cnsh r 

dr cosnr r (2Zar) _ j + 2 (£ a )V 2 

sinh r 

[1 - cos(2Z«r)] = Re[^(l + iZa) + C] . 



(A4) 



In order to calculate the function g(Za) we make a shift p — > p + q + in (|2^|) and use the 
exponential parameterization 



.4' 



7/2 r ds 



(A5) 



Then we have 
g(Za) 



■ irZa/2 r poo 

4 * qT{iZa) J d PP~ 1+2lZ J ds s-^ Za {eM^P ~ qf] - eM*sp 2 ]} ■ (A6) 



Taking the integral over the angles of two-dimensional vector p and making the substitutions 
P ~^ p/y/s, s — > s/q 2 we come to 

poo 



g(Za) 



ic 



, , dp p 2iZa exp[*p 2 ] / dss-"l 2 [e ls U2p^)-l\ ■ (A7) 
2T(iZa) J Jo 1 
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Taking the integrals over s and then over p we finally obtain 



. * „ r 1 -iZa r (12 + iZa) /A , 

g(Za)=Za-± ') ' — 4 . A8 

yv ; r(l + iZa)r(l/2-iZa) v ; 
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